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ABSTRACT 

We present high resolution N-body/SPH simulations of the interacting cluster 
1E0657— 56. The main and the sub-cluster are modeled using extended cuspy ACDM 
dark matter halos and isothermal /^-profiles for the collisional component. The hot 
gas is initially in hydrostatic equilibrium inside the global potential of the clusters. 
We investigate the X-ray morphology and derive the most likely impact parameters, 
mass ratios and initial relative velocities. We find that the observed displacement be- 
tween the X-ray peaks and the associated mass distribution, the morphology of the 
bow shock, the surface brightness and projected temperature profiles across the shock 
discontinuity can be well reproduced by offset 1:6 encounters where the sub-cluster 
has initial velocity (in the rest frame of the main cluster) close to 2 times the virial 
velocity of the main cluster dark matter halo. A model with the same mass ratio and 
lower velocity (1.5 times the main cluster virial velocity) matches quite well most of 
the observations. However, it does not reproduce the morphology of the main cluster 
peak. Dynamical friction strongly affects the kinematics of the sub-cluster so that the 
low velocity bullet is actually bound to the main system at the end of the simulation. 
We find that a relatively high concentration (c=6) of the main cluster dark matter 
halo is necessary in order to prevent the disruption of the associated X-ray peak. For 
a selected sub-sample of runs we perform a detailed three dimensional analysis follow- 
ing the past, present and future evolution of the interacting systems. In particular, 
we investigate the kinematics of the gas and dark matter components as well as the 
changes in the density profiles and the motion of the system in the Lx — T diagram. 

Key words: methods: N-body simulations - galaxies: clusters: individual: 1E0657-56 
- dark matter - X-rays:galaxies:clusters 



1 INTRODUCTION 

The "bullet" cluster 1E0657-56 represents one of the most 
complex and unusual large-scale structures ever observed. 
Located at a redshift z — 0.296 it has the highest X-ray 
luminosity and temperature of all known clusters as a re- 
sult of overheating due to a recent supe rsonic Mach M ~ 3 
(JMarkevitch et al.ll2002l : lMarkevitchll2006. ) central encounter 
of a sub-cluster (the bullet) with its main cluster. The 500 
ks Chandra ACIS-I image of 1E0657-56 (Fig. 1 in Marke- 
vitch 2006) shows two plasma concentrations with the bul- 
let sub-cluster on the right of the image being deformed in 
a classical bow shock on the western side as a result of its 
motion through the hot gas of the main cluster. The anal- 
ysis of the shock structure leads to the conclusion that the 
bullet is now moving aw ay from the main cluster with a ve- 
locity of - 4700 km s"^ (|Markevitchll2006l ) . The line-of-sight 
velocity difference between the two systems is only 600 km 
s~^ suggesting that the encounter is seen nearly in the plane 
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of the sky ([Barrena et al.l [2002) • As the core passage must 
have occurred ~ 0.15 Gyr ago we have the unique oppor- 
tunity to study this interaction in a very special short-lived 
stage, far away from thermal and dynamical equilibrium. 
As a result of the encounter, the collision-dominated hot 
plasma and the collisionless stellar and dark matter com- 
ponents have been separated. The galaxy components of 
both clusters are cl early offset from the associated X-ra y 
emitting cluster gas (JLiang et al.ll2000l : iBarrena et al.ll200^V 
In addition, weak and strong lensing maps (|Clowe et al.l 
12004 I2OO6I : iBradac et al]l2006l ) show that the gravitational 
potential does not trace the distribution of the hot clus- 
ter gas that dominates the baryonic mass but follows ap- 
proximately the galaxy distribution as expected for a col- 
lisionless dark matter component. The likelihood to find 
such a high velocity cluster encounter in a ACDM cos- 
mology has recently been investigated bv lHavashi fc White! 
(|2006 ) using the Millenium Run simulation . According to 
the newest est imates from X-ray and gravitational lens- 
ing results the iHavashi &: Whitd (120061) likelihood becomes 
-- 0.8 X lO"'^ (JFarrar k RoseiJ l2007l ) which means that 
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1E0657-56 represents an extremel y rare system in a ACDM 
universe. Recent numeri cal works (|MilQsavlievic et al.ll2007l : 
ISpr ingel & Farrarl 120071 ) have demonstrated that the rela- 
tive velocity of the dark matter components associated to 
the main and the sub-cluster are not necessarily coincident 
with th e bullet velocity inf e rred f rom the shock analysis. In 
details. iMilosavlievic et al.l (|2007h using a 2-D Eulerian code 
well reproduced the observed increase in temperature across 
the shock front with a dark halo velocity^ 16% lower than 
that of the shock, while ISpringel k, Farrarl (|2007h found even 
a larger difference between the shock velocity (~ 4500 km 
s~^ in their best model) and the spee d of the halo (only 
^ 260 km s~^). Moreover, according to IMilosavlievic et al.l 
(|2007|), due to a drop in ram-pressure after the cores' inter- 
action the gas component of the sub-cluster can eventually 
be larger than that of its dark matter counterpart. 

The simulations of Springel & Farrar (2007) represent 
the most complete three dimensional numerical modeling of 
the 1E0657-56 system so far. Nevertheless they focus prefer- 
entially on the speed of the bullet but fail in reproducing the 
observed displacement of the X-ray peaks which represent 
an important indicator of the nature of the interaction. In 
particular, they do not obtain any displacement in the X-ray 
distribution associated with the main cluster suggesting that 
the baryonic component is suffering too little ram-pressure. 
Moreover, the concentrations used for the main cluster (and 
obtained by modeling the lensi ng data) are much smaller 
than those suggested by ACDM (JMaccio et al.ll2007[ ) for ha- 
los of similar masses. 

The aim of this paper is to investigate the evolution of 
the bullet cluster in details using high resolution SPH simu- 
lations. We quantify the initial conditions that are required 
in order to better reproduce its observed state and predict 
its subsequent evolution. 

Our model allows us to determine in details the spatial, 
thermal and dynamical state of the dark matter and hot gas 
distribution in the bullet cluster. 

The paper is organized as follows. Section 2 describes 
the adopted cluster models and orbital parameters. In Sec- 
tion 3 we perform a projected analysis of our simulations 
comparing it with the latest X-ray and gravitational lensing 
results. In Section 4 we select some significant models and 
study in details the three dimensional kinematics and mor- 
phology of the interacting systems and their past and future 
evolution with time, as well as the motion of the main cluster 
along the Lx — T diagram. 



2 MODELS 

Both the main and the sub- cluster are two components 
spherical systems modeled assuming a cuspy dark matter 
halo and a distribution of hot gas in hydrostatic equilibrium 
within t he global potential of the cluste r. The dark halo has 
a NFW (|Navarro. Frank fc Whit Jll997l ) profile: 
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where pcrit is the critical density of the universe at the time 
of the halo formation, Ts is a scale radius and 6c the char- 
acteristic halo over density. The virial mass Mvir and radius 
rvir are related by: Mvir — Avirpcrit{^7T/3)rlir, where the 



density contrast Ayir is set equal to 200. The concentra- 
tion parameter c = Vyir/rs is assume d to be dependent 
on the halo mass (JMaccio et al.l l2007l l . The velocity dis- 
tribution at a given point in space is approximated by a 
Gaussian, whose velocity dispersion is gi ven by the solu - 
tion of the Jeans equation at this point (|Hernquistlll993l ). 
The distribution of hot gas follow s an isothermal /3- model 
(| Cavalier e fc Fusco-Femianolll976l l of the form: 



p(r) = po(l + ir/r.f) 
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We take the asympt otic slope param eter = 2/3 
(Ijones &: Formanlll984[ ) and Tc — l/2rs (iRicker fc SarazinI 
|200l|). The adopted gas fraction ranges from a minimum 
value of 8 % to 17%. consistent with the recent WMAP re- 
sults (Spe rgel et al.l 12007 ). with an intermediate value of 
12%, comparable with the gas mass fraction provided by 
X-ray observations of galaxy clusters (|Vikhlinin et al.ll2006l : 
iMcCarthvet al.ll2007l ). 

Assuming a spherically symmetric model, the temper- 
ature profile is determined by the condition of hydrostatic 
equilibrium by the cumulativ e total mass distribution and 
the density profile of the gas ([Mastropietro et al.ll2005l ): 
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where M(r) is the total mass within the radius r, rup is the 
proton mass and p the mean molecular weight. We assume 
/i = 0.6 for a gas of primordial composition, which appears 
to be a reasonable approximation since the mean tempera- 
tureof 1E0657-56 is T -- 14 keV according to iMarkevitclJ 
( 120061 ) and cooling is dominated by bremsstrahlung and al- 
most independent of the metallicity. G and ks are the grav- 
itational and Boltzmann constants. 

Masses are assigned to the models acc ording to the 
weak and strong lensing mass reconstruction of lBradac et al.l 
(;2006i V In particular we assume that the inferred mass en- 
closed within t he field of the HST Advanced Camera for 
Surveys (ACS) ([Bradac et al.l 1200 6) is comparable with the 
total projected mass of our simulated system (calculated 
when the two centers of the mass distribution are at a dis- 
tance similar to the observed one) within the same area. 
Since the ACS field represents only the central fraction of 
the area covered by the entire system, this mass constraint 
is strongly influenced by the concentration of the dark mat- 
ter halos (JNusserl 2007). With a cos mologically motivated 
choice of c = 6 ([Maccio et al.l |2007|) for the main cluster 
initial halo model, we can reproduce the lensing mass rea- 
sonably well adopting a main cluster total mass (within the 
virial radius) of ^ 8.34 x lO^^M© (Table [1]), almost a factor 
1.8 smaller than the mass obtained by fltting lensing data 
with extremely low concentrated (c < 2) NFW halos where 
the inner density profile is much fiatter than the one sug- 
gested by ACDM simulations. One of the simulations pre- 
sented in this paper adopts a main halo with a lower (c = 4) 
concentration value. We will see that in this case the X- 
ray intensity peak associated with the main cluster is easily 
destroyed during the interaction. 

We model encounters with mass ratios 1:3, 1:6 and 1:8 
between the sub and the main cluster in order to investigate 
the effects of tidal and ram-pressure stripping, which sig- 
nificantly reduce the mass associated to the sub-cluster and 
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Table 1. Initial conditions of the simulations. For each run dark matter virial mass (M^ir) and concentration c of the main and the 
sub-cluster are indicated, fg, 6, Vi and v are the adopted gas mass fraction, impact parameter and initial velocity of the sub-cluster in 
the system of reference where the main-cluster is at rest and in the center of mass rest-frame, respectively. The last column indicates if 
radiative cooling is included. 



Run 


M^^^(main) 
[1014 Mq] 


c (main) 


M^i^ (bullet) 
[1014 Mo] 


c (bullet) 


/. 


b 

[kpc] 


Vi 

[km/s] 


V 

[km/s] 


Model 


l:6b0 


7.13 


6 


1.14 


8 


0.17 





5000 


4286 


adiabatic 




6v3000b0 


7.13 


6 


1.14 


8 


0.17 





3000 


2571 


adiabatic 




6 


7.13 


6 


1.14 


8 


0.17 


150 


5000 


4286 


adiabatic 




6v3000 


7.13 


6 


1.14 


8 


0.17 


150 


3000 


2571 


adiabatic 




6v2000 


7.13 


6 


1.14 


8 


0.17 


150 


2000 


1714 


adiabatic 




3 


7.13 


6 


2.4 


7 


0.17 


150 


5000 


3750 


adiabatic 




8 


7.13 


6 


0.91 


8 


0.17 


150 


5000 


4445 


adiabatic 




6v3000big 


14.2 


6 


2.4 


7 


0.17 


150 


3000 


2571 


adiabatic 




6c4 


7.13 


4 


1.14 


8 


0.17 


150 


5000 


4286 


adiabatic 




61fg 


7.13 


6 


1.14 


8 


0.12 


150 


5000 


4286 


adiabatic 




31fg 


7.13 


6 


2.4 


7 


0.12 


150 


5000 


3750 


adiabatic 




6c 


7.13 


6 


1.14 


8 


0.12 


150 


5000 


4286 


cooling 




3clfg 


7.13 


6 


2.4 


7 


0.12 


150 


5000 


3750 


cooling 



Table 2. Present time. A is the projected (perpendicular to the plane of the encounter) distance between the peaks of the total mass 
distributions, associated with the two clusters. The third and fourth columns represent the projected offset between each X-ray peak 
and the associated mass density peak. Vgas and v^ark are the sub-cluster gas and dark matter velocity calculated in the center of mass 
system of reference. 



Run 



A offset (bullet) offset (main) 



Vdark 





6b0 


753 


278 


188 


3215 


4715 




6v3000b0 


741 


213 


66 


3131 


3134 




6 


742 


237 


128 


3609 


4756 




6v3000 


729 


185 


172 


2893 


3137 




6v2000 


721 


126 


230 


2849 


2425 




3 


784 


162 


223 


3908 


4076 




8 


737 


228 


117 


3647 


4858 




6v3000big 


725 


139 


92 


3927 


3528 




6c4 


735 


151 


192 


4168 


4799 




61fg 


718 


200 


189 


3811 


4804 




31fg 


779 


197 


242 


3746 


4145 




6c 


736 


234 


127 


3497 


4806 




3clfg 


780 


228 


272 


3595 


4205 



lead to value s closer to the 1:10 ratio inferred from lensi ng 
observations (jClowe et al.ll2004l . l2006l : iBradac et al.ll2006h . 

The main cluster is initially at rest and the sub-cluster 
moves in the x direction with a velocity which ranges from 
2000 to 5000 km s~^. The initial conditions of the different 
runs are summarized in Table [1] The velocities of the sub- 
cluster relative to the center of mass of the system are listed 
in the last but one column of the Table. 

All the simulations were carried out usin g GASOLINE, 
a par allel SPH tree-code with multi-stepping ([Wadslev et al.l 
|200J). Most of the runs are adiabatic, with 7 = 5/3. Radia- 
tive cooling for a primordial mixture of hydrogen and helium 
in collisional equilibrium is implemented in 1:6c and l:31fgc. 
The consequences of assuming a lower concentration for the 
main dark halo are investigated in run l:6c4, where c = 4. 

The main cluster is modeled with 1.8 x 10^ particles, 
10^ SPH and 8 x 10^ collisionless. The sub-cluster, with the 
exception of run 1:8 (where the number of gas particles in 
the bullet is 4 x 10^), has 9 x 10^ particles, 5 x 10^ collisional 



and the remainder dark matter particles. The gravitational 
spline softening is set equal to 5 kpc for the gaseous and 
dark component. 



3 PROJECTED ANALYSIS 

Each numerical work which aims to simulate the bullet clus- 
ter should be able to reproduce simultaneously the main 
features observed in X-ray maps (the bow shock, the rela- 
tive surface brightness of the bullet and the main cluster), 
and the observed surface brightness and temperature pro- 
files across the shock discontinuity. An additional constraint 
is provided by the observed displacement between X-ray and 
lensing maps, which is not negligible in both the main and 
the sub-cluster (^ 110 and ^ 270 kpc respectively, according 
to Clowe et al. 2006). 

A first indication about the validity of a model arises 
from the qualitative comparison of our simulated X-ray sur- 
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Figure 1. Upper panel. 500 ks Chandra image of the system 
with weak lensing k reconstruction shown in green (courtesy of D. 
Clowe). Central and bottom panel. 0.8 — 4 keV surface brightness 
maps of runs l:6vb0 and I:6v3000b0. Logarithmic colour scaling 
is indicated by the key at the bottom of the figure with violet 
corresponding to 10"^^ erg s"-*^ kpc ~^ and white to 2 x lO^-*^ erg 
s"-*^ kpc ~^. White contours trace the total surface mass density 
of the system within 2.3 x 10^ and 2.3 x 10^ M© kpc~2. The box 
size is 1800 kpc. 



face brightness map s with the X - ray 50 ks Chandra ACIS-I 
images provided bv iMarkevitchI (|2006h . 

The impact parameter h is not strictly constrained by 
observations. Nevertheless a head-on merger, with 6 = 0, 
seems to be excluded by comparing deep X-ray observations 
and weak lensing maps. In particular i n the top panel of Fig. 
1 - which is Fig. lb of lClowe etHI (|2006h - the brightest 
gas associated with the main cluster is not located along the 
line which connects the centers of the two total mass distri- 
butions. Moreover, the X-ray emission from the main cluster 
is asymmetric, with a peak in the north of the image and an 
extended tail of less bright material pointing south. These 
features are hardly associable with a zero impact parameter 
interaction, as shown in the middle and bottom panels of 
Fig. [1] which illustrate two simulated 1:6 head-on encoun- 
ters where the sub-cluster moves from the left to the right 
of the image (x axis of the simulation) with decreasing ini- 
tial velocities (the middle panel corresponds to run l:6vb0 
of Table [1] the bottom one to I:6v3000b0). Images are pro- 
jected along an axis perpendicular to the collision plane (the 
encounter is seen face-on) and the selected snapshot is the 
one which most closely matches the observed distance be- 
tween the centers of the total mass distributio ns, associated 
with t he two clusters - about 720 kpc from iBradac et al.l 
((20061) - once the sub-cluster has passed through the core of 
the main system (hereafter the present time). 

In the lower two panels of Fig. [1] colours represent X- 
ray maps in the Chandra energy band (0.8-4 keV) generated 
using the Theoretical Image Processing System (TIPSY), 
which produces projected X-ray surface brightness maps 
with the appropriate variable SPH kernel applied individ- 
ually to the flux represented by each particle. Assuming 
complete ionization and zero metallicity (metal lines are 
expected to provide a significant contribution to emissiv- 
ity only at relatively low temperatures, less than 2 keV ) 
the X-ray luminosity in a given energy band is defined as 



tne A-ray luminosity 
(|Borgani et al.ll2004l ) 



Lx = (m^p) 



Ngas 



mipiA{Ti) 



(4) 



where A{Ti) is the cooUng function in the specific band, T^, 
pi and rrii are temperature, density and mass associated with 
the i — th hot (Ti > 10^ K ) gas particle, respectively, nip is 
the proton mass and p— 0.6 the mean molecular weight. The 
sum runs over all the Ngas particles within an oblong of base 
equal to the pixel size and major axis oriented along the line 
of sight. When the X-ray luminosity of the entire cluster is 
calculated Ngas is the number of hot particles within the 
virial radius rvir- The c ooling function is compu ted using 
a Raymond-Smith code ((Ravmond fc Smithlll977h for a gas 
of primordial composition. The energy band (0.8-4 keV) is 
chosen in such a way to reproduce the 500 ks Chandra ACIS- 
I image of the bullet cluster ()Markevitchll2006l ) . The entire 
energy band, used to calculate bolometric X-ray luminosities 
in the following of the paper, goes from 5 eV to 5 x 10"^ keV. 
The lower two panels of Fig. [1] show that after an en- 
counter with zero impact parameter the displacement of the 
main cluster's X-ray peak is aligned with the x axis. A large 
relative velocity (central panel) induces a significant offset 
between the dark and baryonic component of the main clus- 
ter (see Table [2] for details) but it also leads to substantial 



Simulating the Bullet Cluster 5 



disruption of the main cluster gaseous core. Moreover, the 
displacement of the bullet from its dark halo (278 kpc) is 
much larger than observed. Decreasing the relative velocity 
between the two clusters (bottom panel) two X-ray peaks 
are clearly visible but the displacement of the main cluster 
gas is now negligible due to the lower ram- pressure expe- 
rienced by the main cluster core. ISpringel fc Farrarl (|2007|) 
(Fig. 7 in their article) provide further examples of head-on 
encounters with even lower mass ratios and relative veloc- 
ities {v = 2600 km s""*^ in the center of mass rest frame). 
Even assuming extremely low concentrations (c = 2) for the 
main halo, the authors never reproduce the displacement 
observed in the two systems. Increasing the concentration 
strongly increases the luminosity of the main cluster, which 
appears much brighter than the bullet, contrary to what is 
observed. 

A bow shock is clearly visible on the right of each image. 
The shape of the shock front is only marginally dependent 
on the kinematics of the model while the distance between 
the edge of the bullet (the so called contact discontinuity) 
and the shock front becomes larger for decreasing bullet ve- 
locities. The contact discontinuity itself is much flatter in the 
case of l:6b0 than in the low velocity encounter I:6v3000b0 
and clearly not comparable with observations, which show 
a more narrow structure. In general, a more efficient ram- 
pressure during the phase of core-core interaction is associ- 
ated with a larger op ening an^le of the con tact discontinuity 
at the present time (JQuilis k Moorell200il ). 

The rest of the runs listed in Table [1] have an impact 
parameter b equal to 150 kpc, comparable with the core 
radius Tc of the main cluster gas distribution for most of the 
models. For a plasma distributed according to a /3 profile 
like the one adopted in this paper the assumption 6 = 150 
kpc implies that the maximum external density crossed by 
the core of the sub-cluster is 1.5 times smaller than the one 
it would pass through for 6 = kpc. For the same initial 
velocity, the sub-cluster also sees a weaker potential with 
respect to the case of a head-on encounter and experiences a 
smaller maximum orbital velocity. A choice of a much larger 
value of b would further decrease the mutual ram-pressure 
between the two systems and require much higher relative 
velocities in order to explain the observed offset between the 
gas and dark matter. 

In Fig. [2] and [3] we illustrate the projected X-ray surface 
brightness maps of some interesting models with impact pa- 
rameter b = 150 kpc. The encounters are shown face-on. The 
box size and the surface mass density contours are the same 
as in the last two panels of Fig. [1] In order to underline the 
morphological details of the high emission regions the upper 
limit of the surface brightness scale varies in the different 
images. Individual values are indicated in the captions. 

Among this sub-sample of runs, 1:3 produces the largest 
displacement of the X-ray peak associated with the main 
cluster, but the X-ray map differs from the observations. In 
particular a large strip of strongly emitting gas still connects 
the two X-ray peaks while the morphology of the main clus- 
ter peak is much more elongated than observed. Decreasing 
the mass ratio between the two interacting systems (runs 1:6 
and 1:8) the displacement in the bullet becomes larger than 
that in the main cluster. The run 1:6 is characterized by an 
initial sub-cluster velocity of 5000 km s~^ (as well as runs 1:3 
and 1:8) in the system of reference where the main cluster is 



at rest, that corresponds to a present time velocity of ^ 4300 
km s~^ in the center of mass rest frame. The same model is 
simulated assuming lower relative velocities (l:6v3000 and 
l:6v2000). With decreasing velocities the offset in the sub- 
cluster becomes smaller and the X-ray emission from the 
bullet less pronounced with respect to the bright X-ray emit- 
ting region at the center of the main cluster. At the same 
time the shape of the contact discontinuity changes, getting 
progressively more narrow while the distance between the 
contact discontinuity and the shock front increases as will 
be shown more quantitatively in the next section. The dis- 
placement associated with the main-cluster is determined by 
the distance of closest approach between the centers of the 
two clusters, which becomes smaller - assuming the same 
initial impact parameter b - with decreasing bullet veloci- 
ties. Indeed, the separation between the main cluster X-ray 
emission peak and its dark matter counterpart is maximum 
in the case of the low velocity run l:6v2000. 

A not negligible fraction of the X-ray emission visible 
at the present time near the center of the main-cluster is 
actually associated with hot gas stripped from the exter- 
nal regions of the bullet. Fig. prefers to run l:6v3000. It 
shows the individual distribution of gas originating from 
the main (top) and sub-cluster (bottom) and lying at the 
present time within 1 Mpc from the center of the system. 
Comparing these images with the middle right panel of Fig. 
[2] it appears evident that the bright elongated X-ray feature 
crossing the second innermost isodensity contour is associ- 
ated with the displaced gaseous center of the main cluster 
while the surrounding more diffuse region hosts a significant 
amount of sub- cluster gas. Indeed, as it shown in the upper 
panel of Fig. 2] the motion of the bullet across the inner 
regions of the main system creates a low density "tunnel" 
in the main cluster gas distribution. At the same time the 
sub- cluster looses a large amount of baryonic material dur- 
ing the phase of core-core interaction. This material, which 
fills the tunnel, falls back into the gravitational center of the 
main cluster and resides now at more than 500 kpc distance 
from the X-ray bullet. As we will show in the next section, 
the amount of gas deposited by the sub-cluster in the cen- 
tral regions of the main system increases with decreasing 
relative velocities. This trend explains the relative increase 
in luminosity of the diffuse strongly emitting component if 
compared to the peak associated to the main cluster core 
gas (always in the upper right region with respect to the 
center of the mass density distribution), when we compare 
run 1:6 with the low velocity run l:6v2000 where it becomes 
the primary peak of X-ray emission. 

As shown in the bottom right panel of Fig. [2] a main 
halo with low concentration (c = 4) does not survive a 1:6 
sub-cluster encounter with velocity v = 5000 km s~^ and its 
X-ray peak is destroyed. 

As previously noticed by ISpringel fc Farrarl (|2007|) the 
choice of a lower gas fraction (fg = 0.12 in three of the last 
four runs of Table [1]) does not affect significantly the X-ray 
map morphology although the displacement of the two lumi- 
nosity peaks with respect to their dark matter counterparts 
slightly changes. 

Finally, we tested the consequences caused by includ- 
ing radiative cooling and choosing a larger main halo model 
(Fig. [3|). Cooling makes the contact discontinuity more 
narrow and the amount of diffuse X-ray gas around the 
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Figure 2. 0.8-4 keV surface brightness maps of runs (from the top left to the bottom right) 1:3, 1:8, 1:6, l:6v3000, l:6v2000 and l:6c4. 
Logarithmic colour scaling is indicated by the key at the bottom of the figure, with violet corresponding to 10^^ ergs~^ kpc~^ and white 
to 1.8 X 10^-*^ erg s~^ kpc"-^ in runs 1:8 and l:6c4, to 2.34 x 10^-*^ erg s~^ kpc"-^ in run 1:3 and to 2 x 10^-*^ ergs"-*^ kpc"-^ in the remaining 
cases. Projected isodensity contours of the total mass distribution are shown. Limits are 2.3 x 10"^ and 2.3 x 10^ M© kpc~^. Each box 
size is 1800 kpc. 



Simulating the Bullet Cluster 7 




Figure 3. Same as in Fig. [2]for the cooling run 1:6c (the X-ray upper hmit is 1.8 x lO"^-*^ erg s -"^kpc ^) and run l:6v3000big (3.9 x lO'^^ 
erg s~-'^kpc~^). 



peak associated with the main cluster smaller. The off- 
set in the main and sub-cluster remains unaltered. Run 
l:6v3000big is characterized by a main cluster total mass 
of ^ 1.64 X 10^^ Mq, which is closer to the value adopted 
by ([Springel fc Farrarl l2007l l and predicted by fitting the 
large field weak lensing data with extremely low concen- 
trated (c ~ 2) NEW halos. Although the initial relative 
velocity is only 3000 km s~^, due to the large mass of the 
host halo the present time velocity of the bullet in the center 
of mass rest frame is much higher than in the correspond- 
ing l:6v3000 run. Consequently the offset of the X-ray peak 
associated with the main cluster is less than 100 kpc and 
the amount of gas lost by the bullet in the core of the main 
system closer to that observed in 1:6 than in l:6v3000. 

Summarizing, run l:6v3000, with mass ratio 1:6, initial 
relative velocity v = 3000 km s~^ and present time sub- 
cluster velocity v ^ 3100 km s~^ in the center of mass rest 
frame, best reproduces the main features observed in X- 
ray maps, in particular the peculiar morphology of the X- 
ray emission associated with the main cluster, the relative 
surface brightness between the main and the sub-cluster, the 
shape of the shock front and of the contact discontinuity. 
Although the low velocity run l:6v200 Q leads to a X-ra; 
displacement closer to that observed bv lClowe et al.l (|200' 
both in the main and in the sub-cluster, this model seems to 
be excluded on the basis of a pure morphological comparison 
with the observational data. Indeed, the bullet seems to be 
much less bright than the center of the main cluster. 

Fig. [5] compares the integrated mass profiles of the main 
and the sub- cluster for differen t mass mode l s with the grav- 
itational lensing res ults from ,Bradac et al. ([20061 ) and pre- 
vious simulations by ISpringel fc Farrarl (|2007|) . The profiles 
inferred from observations are obtained by measuring the 
enclosed mass in cylinders centered on the southern cD of 
the main cluster and on the BCG of the sub-cluster. Here we 
calculated the present time projected mass within cylindri- 
cal bins centered on the center of mass of the two clusters. 
The massive cluster run l:6v3000big (green lines), overesti- 



mates the cumulative profile of the main cluster, while an 
initial main cluster mass of ^ 8.3 x 10^^ Mq (1:6v3000, in 
red) fits the observational p oints better than the larg e-mass- 
low-concentration model of ISpringel fc Farrarl (|2007|) in the 
central regions and underestimates the projected mass by 
less than 20% between 250 and 450 kpc from the center. Blue 
lines refer to run l:31fg where the main halo has the same 
dark matter mass than l:6v3000 but lower baryonic fraction. 
It is interesting to notice that using almost the same initial 
sub-cluster mass (the sub-cluster model in run l:6v3000big 
and l:31fg only differs in the baryonic fraction) we can nearly 
reproduce the observed projected mass in the bullet only 
in the case of a very massive host halo (green dashed line, 
l:6v3000big run) which overpredicts the main cluster profile 
(green solid line). On the other hand a smaller main halo 

- which better fits the observational data (blue solid line) 

- associated with a similar sub-cluster mass (l:31fg) under- 
predicts the lensing results (blue dashed line). Our favourite 
run l:6v3000 gives even lower projected sub-cluster mass 
estimates (red dashed line) and lies close to the model pro- 
posed by Springel & Farrar (2007) (their sub-cluster mass 
is similar to the one used in our 1:6 models). These results 
seem to indicate that the mass of the sub-cluster is a signif- 
icant part of the main cluster mass, with initial mass ratio 
between the s ub and the main cluster much larger than 1:3. 
NusseJ (|2007|) successfully reproduces the projected mass 
associated with the sub-cluster only assuming very massive 
halos and high mass ratios. In particular the smallest mass 
ratio allowed by their models is 1:2.7, with a main cluster 
and bullet mass of 3.2 x lO^^M© and 1.2 x lO^^M©, respec- 
tively. Such high masses are however clearly incompatible 
with the values inferred b y galaxy kinematics. In particular 
the virial mass derived bv Barre na et al.l (|2002n for the main 
cluster is 1.24 x 10^^ (comparable with the value obtained 
by Girardi & Mezzetti 2001). They also measure the main 
cluster total R-band luminosity and f ind Lr — 10^^ Lc?). Us- 
ing Lr to estimate the cluster mass (JMiller et al.ll2005l ) we 
obtain M200 '^ lO^^M^. The bullet mass determination of 
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Figure 4. Run l:6v3000. Gas originating from the main (top 
panel) and the sub-cluster (bottom) is projected individually 
along the z-axis perpendicular to the plane of the encounter. Vi- 
olet corresponds to a surface density of 2.3 x 10^ M© kpc~^ and 
white to 2.3 x 10^ M© kpc~-^. Projected isodensity contours of the 
total mass distribution are drawn on top of the image. The box 
size is 1 Mpc. 



iBarrena et alj (|2002l l is less reliable since it is based on only 
seven galaxies and on the assumption of equilibrium while 
the sub-cluster seems to be tidally perturbed beyond ~ 200 
kpc from its center. They find M200 = 1.2 x 10^^ M©. From 
the total R-band luminosity {Lr — 0.2 x lO^^L©) we get a 
value almost five times l arger (M200 ^ 5 x 10^^ Mq). The 
bullet mass proposed by INusserl ( 20071 ) is even one order 
of magnitude larger than that obtained from weak lensing 
analysis (1.5 x 10^ ^0) of a new larger lensing field which 
covers most of the area occupied by the system (Clowe et al. 
in prep.). With such low mass ratios we cannot reproduce 
the morphology of the X-ray maps (upper left panel of Fig. 
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Figure 5. Cumulative pro j ect ed mass profi les . Points represent 
the observational results (Brad ac et al.|[2006l ) for the main (dots) 
and the sub-clust e r (sta rs). Black curves refer to the model of 
ISpringel fc Farran (|2007l ) (solid and dashed black line indicates 
the main and sub-cluster mass, respectively). Coloured lines refer 
to the models in this work, always with solid curves indicating the 
main cluster and dashed ones the bullet. In details the red lines 
are associated to our favourite model l:6v3000, blue to l:31fg and 
green to l:6v3000big. 



[2|) since the luminosity peak associated with the main cluster 
core is almost completely destroyed by the interaction. 

Actually any attempt to derive the mass of the main 
and sub- cluster simply by projecting the mass of two iso- 
lated NFW halos and fitting t he observed lensing data - 
like the analysis of lNusserl (|2007|) - is going to provide wrong 
results. Indeed, as we will show in Section 4, the dark mat- 
ter profiles of the main and sub-cluster change significantly 
during the central phases of the interaction, with a density 
increase in the inner regions and a decrement beyond 0.3 
Tvir- This effect, which characterizes both the main and the 
sub-cluster, is not related to stripping processes but can be 
explained by the fact that the two halos are strongly per- 
turbed and cannot be idealized as the sum of two isolated 
equilibrium models. Fitting the projected cumulative mass 
around the centers of the two interacting clusters with NFW 
models we would overestimate either the virial mass or the 
concentration of t he original clusters. 

According to iBarrena et al.l (|2002l l the line of sight ve- 
locity difference between the main and the sub-cluster is 
relatively small, about 600 km s~^, which implies that the 
encounter must be occurring nearly in the plane of the sky. 
In Fig. [6] we show the X-ray surface brightness map of run 
l:6v3000 when the orbital plane is viewed inclined by 45 
degrees left (top panel) and down (bottom panel). In the 
upper image the bullet (which is closer to the observer) ap- 
pears to be much brighter than the main cluster while the 
contact discontinuity and the shock front are much rounder 
than in the case where the encounter is seen perpendicu- 
lar to the orbital plane. Such a large inclination has to be 
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excluded by a comparison with observations. Moreover, the 
Une of sight velocity difference between the two clusters (we 
assume the line of sight kinematics of the dark matter to 
be coincident with that of the stellar component observed 
by Barrena et al. 2002) is much larger than the observed 
one (~ 2700 km s""*^) due to the fact that the main com- 
ponent (~ 3200 km s""*^ in the center of mass rest frame) 
of the sub-cluster velocity is oriented in the direction of the 
positive X axis in the original face-on encounter. A velocity 
difference comparable to the observed one is reached for a 
right to left inclination of ^ 10 degrees. On the opposite, 
a top-down inclination of the plane of the encounter with 
respect to the Une of sight (bottom panel) does not produce 
drastic changes in the morphology of the bow shock and of 
the edge of the bullet. Indeed the main and the sub-cluster 
still lie close to the plane of the sky, but the relative surface 
brightness between the dense and the diffuse strongly emit- 
ting components at the center of the main system changes. 
In particular the peak of X-ray emission due to the pres- 
ence of hot gas deposited by the sub-cluster within the core 
of the main cluster is now comparable with the main clus- 
ter X-ray peak. The difference in line of sight velocity in 
this case is only ^100 km s~^ and would probably not be 
distinguishable with respect to a pure face-on encounter. 

Fig. [71 represents the X-ray brightness profile measured 
in a narrow slit (with thickness 20 kpc) parallel to the x- 
axis across the shock front. The bullet is located at x = 
with X increasing towards the pre-shock region. The general 
trend does not depend on the model a nd is comparable with 
the profile suggested by observations (|Markevitchll2006h : all 
the simulations show an inner bump associated with the 
bullet and an outer one (between x = 20 and x = 100 kpc) 
caused by the shock. The abrupt jump at x '^ 100 kpc is the 
shock front while the pre-shock region is well fitted by a two 
dimensional /3-profile Sxix) = Sxo[l + (x/xc)^]"^^+^/^ In 
the case of run 1:6 Xc = 150 kpc, (3 = 2/3 and Sxo = 8 x 
10^° erg s~^ kpc~^. A different orientation (45°, bottom- up 
inclination of the plane of the encounter) of the line of sight 
does not affect substantially the surface brightness profile. 
If cooUng is activated the surface brightness of the pre and 
post-shock region is smaller, but the thickness of the shock- 
front and the jump in surface brightness are similar. 

In order to calculate projected temperatures we need 
to define a weighting function. The emission weighted tem- 
perature Tew was originally introduced to provide a bet- 
ter comparison between simulations and observations with 
respect to a simple mass weighted temperature definition 
and has been commo nly used in the analysis of simulations 
(Borgani et alJ|200J). It assumes a weighting function pro- 
portional to the emissivity of each hot gas particle and is 
defined as 



T — 



Y,i''""'mipiA{T^)T^ 



(5) 



J2i <"''mipiA{Ti) 

Recently iMazzotta et al.l (|200^) have demonstrated that 
for clusters with a complex thermal structure the emis- 
sion weighted temperature always overestimates the spectro- 
scopic temperature obtained from X-ray observations due to 
fact that the source is not a single temperature plasma. For 
clusters with a temperature T > 3 keV the discrepancy is 
~ 20 — 30%. This difference becomes particularly large in the 
presence of strong temperature gradients like shocks which 




Figure 6. X-ray surface brightness maps for different inclinations 
of run l:6v3000. In the top panel the orbital plane is left rotated 
by 45 degrees with respect to the projection plane while in the 
panel on the bottom the encounter is seen with a top-down in- 
clination of 45 degrees. The colour scaling is the same as in the 
corresponding face-on map (Fig. [2] middle right panel). The box 
size is 1.8 Mpc. 



appear to be much weaker in observations than what is pre- 
dicted by emission-weighted temperature maps (Mathiesen 
k E vrard 2001, Gar d ini et al. 2004, Rasia et al. 2005). 

iMazzotta et al.l (| 20041 ) proposed a new definition of 
temperature, the spectroscopic-like temperature Tsi: 



1/2 






(6) 



When applied to clusters hotter than 2-3 keV this equation, 
with a = 0.75, gives a good approximation (within few per- 
cent) of the spectroscopic temperature obtained from data 
analysis of Chandra. 

The bullet cluster 1E0657-56 is characterized by 
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Figure 8. Spectroscopic-like temperature profiles measured in a narrow slit (20 kpc) across the shock. The bullet is located at x = kpc 
with X increasing towards the pre-shock region. Upper left panel: 1:6 runs with different relative velocities are compared. Upper right: 
different mass ratios. Bottom left: comparison between adiabatic and cooling 1:6 and 1:3 encounters. Bottom right: effects of inclination, 
lower gas fraction and lower concentration of the main halo. 



the highest luminosity, temperat ure and the stro ngest 
bow shock of all known clusters ([Markevitchl 120061) and 
seems indeed to be the ideal candidate for adopting the 
spectroscopic-like temperature definition. Therefore in the 
remainder of this paper Tsi will be analyzed and Tew only 
indicated for a comparison. Fig. [8] illustrates the projected 
spectroscopic-like temperature profiles across the shock for 
the different runs of Table [1] All the values refer to the 
present time. In particular, the upper left panel shows the 
temperature jumps associated with different relative veloc- 
ities of the two clusters. Decreasing the initial relative ve- 
locity from 5000 to 2000 km s~^ (in the system of reference 
where the main cluster is at rest) reduces the temperature 



peak by ~ 7 keV while the peak itself becomes broader since 
the thickness of the shock increases by almost a factor of two 
due to the lower pressure exercised by the pre-shock gas. 
Both the l:6v3000 and l:6v2000 models seem to fit quite 
nicely the observed height (^ 27 — 30 keV) and th ickness 
(150-200 kpc) of the shock front (|Markevitchl l2006l ) while 
the 1:6 run produces a peak which is too narrow (~ 100 kpc) 
and pronounced (^ 35 keV). The upper right panel refers to 
different mass ratios. If the encounter is characterized by the 
same initial relative velocity and gas fraction the strongest 
shock is associated with the most massive sub-cluster. A 
1:3 adiabatic encounter produces a ^ 45 keV temperature 
peak with thickness ^ 150 kpc while the 1:8 temperature 
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Figure 7. Runs 1:6 and 1:6c: X-ray surface brightness profiles 
across the shock discontinuity. The bullet is located at x = kpc 
and the shock front at ~ 100 kpc. 



Figure 9. Run l:6v3000: spectroscopic-like (Tsi), emission 
weighted (Tew) and true (Tghock) temperature profiles across the 
shock. 



profile is not substantially different from the one of the 1:6 
run. In the same plot we also show the shock created in the 
massive run l:6v3000big. Clearly the maximum temperature 
is much higher than the observed one. Including radiative 
coohng (bottom left panel) has the effect of reducing the 
peak in temperature but does not influence the thickness of 
the shock region. The adoption of a simple cooling model 
is actually questionable. Indeed, although the estimate of 
the temperature jump is reasonable, the entire temperature 
profile drops by 5 keV due to the fact that once cooling is ac- 
tivated the main cluster gas component becomes thermally 
unstable in the early phases of the interaction and overcools 
in the central regions. If cooling is important, models with 
higher mass ratios and relative velocities (like the l:31fg run 
in plot) have still to be taken in account and can not be 
excluded as the high temperature peaks could actually cool 
significantly. The choice of a different line of sight (bottom 
right plot) does not affect significantly the temperature pro- 
file along the shock. Even decreasing the baryonic fraction 
in the clusters and assuming a much less concentrated main 
halo, the height and thickness of the temperature peak do 
not change. 

In Fig.[9]we compare the spectroscopic-like temperature 
profile across the shock region with the emission weighted 
one for our favourite run l:6v3000. While the projected tem- 
perature profile calculated according to the two definitions 
is similar (Tew is only slightly higher than Tgi) in the pre 
and post-shock regions, the emission weighted temperature 
Tew in the region ^ x ^ 150 kpc is ^ 20% higher than 
Tsi. Tshock represents the actual temperature along thex-axis 
through the shock. Indeed the blue curve in Fig. [9] gives the 
exact temperature jump across the shock (the "true tem- 
perature" of the shock) which is characterized by an even 
higher peak with respect to the projected ones. This de- 
projected temperature p rofile is actually the one shown by 
ISpringel fc Farrarl (|2007h in their Fig. 9 where they compare 



their model with observ ations which on t he other hand refer 
to projected quantitie s (IMarkevitchlbOOa ) . As t he calculated 
local temperatures of Sprin gel fc Farrarl (|2007|) fit the pro- 
jected observed temperature profile very well we conclude 
that their projected temperatur e profiles are actua lly incon- 
sistent with the observations by JMarkevitchl (|2006l l . 

In Fig.[Tni[II][l2]we show for our selected adiabatic run 
l:6v3000 the evolution with time of the 0.8—4 keV X-ray sur- 
face brightness, spectroscopic-like temperature and gas sur- 
face density during the central phases of the interaction. All 
the quantities are projected along an axis perpendicular to 
the plane of the collision. On top of the color maps we draw 
in black the projected isodensity contours of the total mass 
distribution, which is dominated by the dark matter compo- 
nent. Time increases from the left to the right and from the 
top to the bottom. The sequence of six panels covers an in- 
terval in time of 400 Myr, between 1.3 and 1.7 Gyr from the 
beginning of the simulation. The snapshot corresponding to 
the present time is the one on the bottom left. The bullet 
approaches the main cluster from the left, with an initial 
velocity of 3000 km s~^ in the main clust er rest-frame. The 
shock front has an arc- like shape (Ricke r fc Sarazinll200lh 
and becomes progressively more asymmetric as the bullet 
moves closer to the core of the main cluster. When the cen- 
ters of the two clusters are less than 250 kpc separated, 
ram-pressure becomes effective in producing a displacement 
(visible in the projected density maps as well as in surface 
brightness) between the gaseous core of the bullet and the 
peak of its associated mass distribution. The offset in the 
main cluster is evident only when the bow shock passes 
through its core. In the central panels of Fig. [10] and [12] 
the X-ray luminosity and surface density saturate in order 
to distinguish features in the maps at later times. Neverthe- 
less, it is visible how the main cluster core gets compressed 
and displaced from the center of the potential towards the 
top right of the image and appears in the X-ray maps (last 
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Figure 10. Time evolution of run l:6v3000. X-ray (0.8-4 keV) surface brightness maps. Logarithmic colour scaling is indicated by the 
key to the bottom of the figure, with violet corresponding to 10^^ ergs"-*^ kpc~^ and white to 2.5 x 10^^ erg s~-'^kpc~^ . Projected 
isodensity contours of the total mass distribution are drawn on top of temperature maps. 
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Figure 11. Time evolution of run l:6v3000. Spectroscopic-like X-ray temperature maps. Logarithmic colour scaling is indicated by the 
key to the bottom of the figure, with violet corresponding to 0.9 keV and white to 86 keV. Projected isodensity contours of the total 
mass distribution are drawn on top of temperature maps. 
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Figure 12. Time evolution of run l:6v3000. The colormap represents the gas surface density projected along an axis perpendicular to 
the plane of the encounter. Violet corresponds to a surface density of 2.3 x 10^ Mq kpc~^ and white to 2.3 x 10^ Mq kpc~^. Projected 
isodensity contours of the total mass distribution are drawn on top of the image. 
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two panels) as an elongated structure, characterized by high 
surface brightness. 

Temperature maps better describe the evolution of the 
shock region, which gets compressed and hotter during the 
core passage while in a later phase it cools down and be- 
comes thicker due to lower pressure of the pre-shock gas. 
The bullet itself expands as it leav es the centra l regio ns of 
the main cluster. As observed by Markevitch (2006) and 
previously noticed bv ISpringel fc F arrar (2007), despite of 
its strong X-ray emission (it is the brightest feature in the 
post core-core interaction X-ray maps) the bullet remains 
relatively cold. Even if radiative cooling is not activated, 
the core of the sub-cluster is heated only to a maximum 
temperature of 10^ K, while the shock front is much hot- 
ter (5 X 10^ K). The projected temperature associated with 
the bullet is higher when it passes through the core of the 
main system (central panels of Fig. Ill|) due to the line of 
sight overlap with the hot gas from the main cluster and the 
shock heated material stripped from the bullet itself which 
surrounds it. As soon as the sub-cluster moves out into the 
cool external regions of the main cluster and loses part of 
the hot envelope of stripped gas, its projected temperature 
decreases to values comparable to the observed ones. An- 
other peculiar feature in the temperature maps is the high 
temperature region next to the innermost total density con- 
tour of the main cluster and visible in the middle right and 
bottom left panels of Fig. [11] . This area could be associated 
with the southeas t ern h i^h temperature region observed by 
iMarkevitch et al.1 (|2002l ) (regions /, i and / in their Fig. 2) 
in the main cluster X-ray map and assumed to be coinci- 
dent with the main merger site. In our simulations the high 
temperature region is filled with hot gas stripped from the 
external regions of the sub-cluster and deposited within the 
core of the main system (compare with Fig. |4]). This high 
temperature material, combined with the main cluster gas 
which lies in the same projected region, produces the diffuse 
X-ray high emission feature visible at the present time in 
the main cluster below the primary peak (bottom left panel 
of Fig. [To] or middle right panel of Fig. [2] for a better colour 
contrast). The eastern side of the high temperature region is 
less bright in X-ray emission and lies beyond the luminosity 
peak associated with the main cluster. As the sub-cluster 
moves to larger radii, the high temperature region expands 
and cools. Shortly (^ 70 Myr) after the present time this 
gas is not anymore clearly distinguishable in the X-ray and 
density maps. 



4 3D ANALYSIS 

In this Section we will focus on a sample of runs from Table 
[1]- namely 1:6, l:6v3000, l:6v2000 and l:3clfg - and in- 
vestigate in details their three-dimensional characteristics. 
Three of these simulations are adiabatic models character- 
ized by different initial relative velocities and will permit us 
to study the effects of the sub-cluster speed on the present 
and future state of the encounter. The remaining run l:3clfg 
is the one which better reproduces - together with l:6v3000 
and l:6v2000 - the observed jump in temperature across the 
shock front. 

Fig.[T3]illustrates the behavior of hydrodynamical quan- 
tities across the shock discontinuity for a snapshot which 



corresponds to the present time. The horizontal axis is cen- 
tered on the gaseous bullet and oriented perpendicular to 
the bow shock nose. The bow shock location is indicated 
by a vertical line while the edge of the bullet corresponds to 
the peak in density. The sub-cluster is moving outward from 
the main cluster core towards positive values of x. All the 
physical quantities are calculated and mass averaged on a 
one dimensional grid where each grid element has a volume 
of 15^ kpc^. V is the component of the velocity perpendicu- 
lar to the shock front and is calculated with respect to the 
system center of mass. The Mach number A^ is determined 
from the temperature jump - which shows a better defined 
discontinuity compared with the density jump - using the 
Rankine-Hugoniot condi tions. We find val ues in good agree- 
ment with observations (JMarkevitchll2006l l for l:6v3000 and 
the cooling run l:3clfg (both with A4 ^ 3). I:6v2000 is char- 
acterized by a slightly lower value of M {^ 2.9) while the 
high velocity run 1:6 has a stronger shock, with M = 3.2. 

As seen in the previous Section, the pre-shock temper- 
ature (~ 11-12 keV for the adiabatic runs) is slightly higher 
than the one obtained by projection along the line of sight 
(Fig. [9] ). Assuming T=ll keV and TW = 3 we predict a 
pre-shock sound speed cs = 1700 km s~^ and a shock ve- 
locity vs = Mcs ^ 5100 km s~^. This value reduces to the 
observed shock velocity vs ^^ 4700 km s~^ if we use the pro- 
jected average pre-shock temperature (T 
bv Markevitch (2006). 

As previous ly noticed by 



("Springel k Fa rraj l2007i ) the velocity 
the last panel of each plot is much 



9 keV) adopted 



other 



authors 
jump shown in 
smaller than the 



theoretically inferred shock velocity. Actually the pre-shock 
gas is not at rest but shows a negative velocity along the 
X-axis which however can only be partially explained by the 
fact that the center of mass of the system is moving in the 
positive direction of the x-axis following the bullet. Indeed 
the upstream velocity maintains a negative sign even with 
respect to the rest frame of the parent cluster, indicating a 
pr e-shock infall towar d s the bullet. This effect is explained 
in Spring el fc Farr ar (2007j by studying the dynamical 
evolution of the system's global potential, which becomes 
deeper after the core- core interaction and induces infall of 
material from the region ahead of the shock. Nevertheless 
the infall velocity which characterizes our models is a sig- 
nificantly smaller than the sub-cluster velocity in contrast 
to the values found by S pringel fc Farrar (2007) and ranges 
only between 500 and 900 km s~^ in the different runs. 

Figs. [14] and [15] illustrate the characteristic velocities 
of the sub-cluster in the orbital plane for the three runs 
l:6v3000, 1:6 and l:6v2000. All the velocities are calculated 
in the center of mass rest frame and the time correspond- 
ing to the present position is indicated by a vertical line. 
The velocity of the dark matter component peaks at the 
moment of closest approach between the two cores and then 
decreases faster than for a ballistic orbit as a result of dy- 
namical friction. The escape velocity at a given sub-cluster 
position is calculated assuming a spherical unperturbed host 
potential and is indicated by a black dotted curve. All the 
runs have initially unbound sub-clusters. Due to the effects 
of dynamical friction after the phase of core-core interac- 
tion the l:6v2000 sub-cluster is actually bound to the main 
system, while sub-clusters with initial velocities v = 3000 
km s~^ and 5000 km s~^ have velocities slightly or much 



16 Mastropietro & Burkert 





60 


> 




CD 


40 


^ 




^ 


20 


CO 





1 
O 


10« 


a 




^ 






10^ 



104 
^ 4000 

!J 2000 

^ 

-2000 






60 


1 1 




> 




(D 


40 


^ 








^ 


20 


CO 





1 
O 


106 


a 




^ 






105 






CL 






lO'' 


1 


4000 




2000 









> 




- 


-2000 




-200 200 

x[kpc] 



400 



-200 200 

x[kpc] 



400 




-200 200 

x[kpc] 



400 



-200 200 

x[kpc] 



400 



Figure 13. Gas temperature, density and one dimensional velocity profiles across the shock discontinuity. The bullet is located at x = 
while the dashed vertical lines indicates the position of the bow shock. 



larger than the escape velocity from the main cluster. The 
gaseous bullet initially follows its dark matter counterpart 
but before the point of closest approach it is slowed down by 
ram-pressure. At the same time the morphology of the sub- 
cluster gas distribution changes. The contact discontinuity 
assumes an arc-like shape (partially reducing the effect of 
ram-pressure) and the bow shock forms. It is interesting to 
notice that after the point of closest approach the evolution 
of the relative velocity between gas and dark matter in the 
sub-cluster strongly depends on the intensity of the ram- 
pressure force. In particular for relatively low ram-pressure 
values (l:6v2000 and l:6v3000) the gaseous bullet is acceler- 
ated towards its dark matter counterpart as soon as it leaves 
the core of the host cluster where it experienced the largest 
external densities and ram-pressure. As a result the grav- 



itational acceleration relative velocity between the gaseous 
and dark component of the sub- cluster (whose two compo- 
nents in the orbital plane are represented for l:6v3000 by 
the blue solid and dashed curves in Fig. [T^ is larger than 
zero. At the time corresponding to the present position the 
two velocities look comparable in the case of l:6v3000 while 
in l:6v2000, where the ram-pressure acting on the bullet is 
lower, the acceleration starts earlier and the gaseous bullet 
is already ^-^ 500 km s~^ faster than its dark counterpart. 
For larger impact velocities - as in the case of 1:6 - ram- 
pressure is effective in slowing down the gaseous bullet even 
at large distances from the center of the main-cluster. The 
velocity of the gaseous bullet is therefore smaller than that 
of its dark counterpart until when the sub-cluster is well 
outside of the virial radius of the main system. 
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Figure 14. Run l:6v3000. Characteristic bullet velocities plotted 
as a function of time. The black and red solid curves represent the 
velocity of the dark matter and gaseous component of the sub- 
cluster in the plane of sky, while the black dotted line indicates 
the escape velocity from the main cluster. The blue curves show 
the relative velocity of the two components of the sub-cluster 
in the two directions perpendicular to the line of sight. The red 
solid and dashed curves show the velocity of the edge of the bullet 
in the center of mass and pre-shock gas rest-frame, respectively. 
Finally, the green curves represent the shock velocity obtained by 
differentiating the shock position. The vertical line indicates the 
present time. More details are given in the text. 
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Figure 15. Run 1:6 (top panel) and l:6v2000 (bottom panel). 
Same as in Fig. [14] 



Figure 16. Bolometric X-ray luminosity L^jol as function of time 
for four selected runs. Luminosity is scaled to its initial value and 
calculated for the entire simulated volume. 



The green solid curve in each plot represents the ve- 
locity of the front shock obtained by differentiating the po- 
sitions of the shock front at increasing times. The shock 
velocity rapidly increases after the core-core interaction and 
at the present time it is always larger than the velocity of 
the gaseous bullet. In order to calculate the bullet and shock 
velocities in the system of reference of the pre-shock gas (red 
and green dashed curve, respectively) the infall velocity of 
the upstream gas is calculated at different times before and 
after the present one. 

The shock velocity vs at the present position is 4100, 
4500 and 5100 km s"^ in the case of l:6v2000, l:6v3000 
and 1:6, respectively, consistently with the shock velocity 
Vs = Mcs inferred from the Rankine-Hugoniot jump 
conditions and with the value provided by observations. 
Only l:6v2000 shows a two-dimensional shock velocity 
well below the o b serva tional uncertainties. As found by 
[Springel & Farrarl (|2007l ), after the point of closest interac- 
tion the shock velocity is always larger than the velocity 
of the sub- cluster mass centroid, but the amount of the 
difference strongly depends on the model. In particular, in 
the case of the two low velocity runs l:6v2000 and l:6v3000 
Vs is ~ 65% and ^ 40% larger than the velocity of the dark 
matter component, while in 1:6 the difference is almost 
negligible (only 6%). 

The cluster collision produces a drastic increase in lumi- 
nosity and temperature. Figs. ll6landll7lshow the bolometric 
X-ray luminosity Lxboi and average spectroscopic-like tem- 
perature Tsi as functions of time. Both quantities are calcu- 
lated for the entire simulated box and scaled to their initial 
values. The first phase of the interaction - which involves 
only the external low density regions of the two clusters - is 
characterized by an identical slow increment of luminosity 
and temperature in all the adiabatic runs. The jump is as- 
sociated with the phase of core-core interaction and peaks 



18 Mastropietro & Burkert 



1 — \ — \ — \ — \ — I — I — \ — \ — \ — \ — I — I — \ — \ — \ — \ — I — r 

l:6v3000 

l:6v2000 

1:6 

l:3clfg 




J_ 



J_ 



50 



2 3 

Time [Gyr] 






10 



J 



1 r^rr 



"1 \ r 

.0 Gyr 

□ present 

^ Markevitch 2006 




l:6v3000 
l:6v2000 
1:6 
l:3clfg 



6 7 8 9 10 
T.1 [keV] 



20 



Figure 17. Average spectroscopic-like temperature Tsi as func- 
tion of time for four selected runs. Temperature is scaled to its 
initial value and averaged for the entire simulated volume. 



right after the time of closest approach. The present time 
is indicated with an empty squared and sits on the down- 
turning curve. The high velocity run 1:6 is associated with 
the largest increase in temperature {Tspec/Tspec{^) ^ 3.5) 
and with the smallest jump in luminosity {Lhoi / Lhoi{^) '^4). 
For decreasing sub-cluster velocities the peak in temperature 
becomes slightly less pronounced while the luminosity jump 
rises by a factor of 1.5. The loss of baryonic material from 
the sub- cluster within the high density core of the main sys- 
tem is indeed larger for low velocity encounters and leads to 
higher luminosities even if the increment in temperature is 
smaller with respect to the high velocity runs. Excluding the 
bound run l:6v2000, the amount of gas stripped from the 
bullet and lying within the virial radius of the main cluster is 
1.6 X 10^^ M© in 1:6 and 1.8 x 10^^ M© in l:6v3000, which in 
both cases corresponds to almost ^ 60% of the initial bary- 
onic content of the sub-cluster. The difference becomes more 
pronounced if we consider only the core {r < Vs) of the main 
cluster, where the mass of sub-cluster gas is 1.4 x 10^^ M© in 
the l:6v3000 run and one order of magnitude less in the case 
of 1:6. At 4 Gyrs after the beginning of the simulation the 
total luminosity of the system is similar to the initial one in 
the case of the bound system l:6v2000 where the center of 
mass of the bullet does not move out to distances beyond 
the virial radius of the main system. The luminosity drops 
to 50% or even less of the initial luminosity for the runs 
l:6v3000 and 1:6 respectively. This decrease in luminosity is 
motivated by the fact that at the final stages of the simu- 
lations a significant fraction of the gas is unbound and very 
extended. Due to its low density it does not contribute to 
the luminosity of the system despite the high temperature. 
In particular, for the same intracluster distance, the high ve- 
locity encounter 1:6 is associated to the highest fraction of 
unbound material, as will be shown later in this Section. On 
the other hand, the final temperature of the system is higher 
than the one associated to the initially isolated clusters and 



Figure 18. Evolution of the interacting system along the Lx —T 
diagram. Here Lxbol cind Tgi are the bolometric X-ray luminosity 
and the spectroscopic-like temperature within the virial radius of 
the main system. The dotted line represents the Lx — T relation 
by Markevitch (1998) for local clusters. The initial and present 
time of each simulation, as well as the observed position of the 
bullet cluster on the Lx — T relation (|Markevitchl l2006l ) are in- 
dicated. 



converges to a value Tsi/Tsi{0) ~ 1.5 almost independently 
of the sub- cluster velocity. The cooling run l:3clfg shows a 
somehow different behaviour: both, luminosity and temper- 
ature profiles have an initial decrement due to the cooling of 
the central regions of the two approaching clusters. Already 
during the early phases of the interaction the cooling run 
moves out of the equilibrium. The peaks in luminosity and 
temperature are much smaller than the corresponding adia- 
batic ones (not illustrated). The final luminosity approaches 

zero^ 

iMarkevitchI (|2006h found that the bullet cl uster Ues ex- 
actly on the Lx — T relation for nearby clusters (JMarkevitchl 
|l998) but its temperature is much higher than the one ex- 
pected according to weak lensing mass estimates. Fig. [18] il- 
lustrates the drift of the simulated systems along the Lx — T 
diagram. The position of 1E0657-560 is indicated by a star. 
Lxhoi and Tsi are calculated in a cylindrical region centered 
on the center of mass of the main cluster and radius equal 
to its initial virial radius. The luminosity of each model is 
normalized in such a way that the initial main cluster lies 
on the Lx — T relation for local clusters despite of the dif- 
ferent initial gas fractions. The starting time of the simula- 
tions is indicated with a black solid square. All the adiabatic 
runs present a similar evolution and move roug hly parallel to 
the L x — T relation, as previously noticed bv lRowlev et al.l 
(| 20041 ) for major mergers in cosmological simulations. Dur- 
ing the early stages of the encounter the cluster moves along 
a curve which is flatter than the observed Lx — T relation: 
the compression of the low density gas at the outskirts of 
the cluster produces an increase in temperature which is 
only marginally accompanied by a luminosity growth. The 
time when the core of the sub-cluster enters the virial ra- 
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Figure 19. Dark matter (upper panels) and gas (bottom) density profiles of the main (left panels) and sub-cluster (right). Initial values 
and profiles at the final time (see text) are shown. Blue curves in the bottom left image refer to gas stripped from the sub-cluster and 
lying in the potential of the main system. Radius is scaled to the virial radius r^ir of the dark matter distribution. 



dius of the main system represents an inversion point in 
the diagram: despite the formation of a bow-shock the tem- 
perature decreases due to the expansion of the main-cluster 
gas and the presence of low temperature baryons belong- 
ing to the sub-cluster within the virial radius of the main 
system. At the same time the luminosity rises as a result of 
the shock and the cluster moves perpendicularly towards the 
Lx — T relation. This phase is actually quite short (on aver- 
age ^0.12 Gyr) but it characterizes all the adiabatic runs. 
Later on the cluster moves almost parallel to the Lx — T 
relation toward larger values of Tgi , with the peak in temper- 
ature being reached at the point of closest approach. Most 
of the runs show a small delay (^ 40 Myr) between the time 
characterized by the highest temperature and the time with 
highest luminosity with the curve making a knot in the di- 
agram. The branch of the curve associated with the post 
core-core interaction is parallel to the increasing one but 
shifted to smaller luminosities: during the strongest phase 



of the interaction some of the hot material is lost beyond 
the virial radius of the main cluster and indeed the largest 
shift in luminosity is observable in the high velocity run 1:6. 
Both, luminosity and temperature now decrease until they 
reach a second inversion point in the curve (in the case of 
1:6 it is only a change in slope), associated with the egress 
of the bullet from the virial radius of the main system. The 
cooling run l:3clfg is characterized by a first decrease in tem- 
perature which corresponds to the initial phase of thermal 
instability and central cooling of the main cluster. Later on 
it moves in the Lx — T diagram similarly to the adiabatic 
runs although the peaks in luminosity and temperature are 
much less prominent. Nevertheless in a pure cooling model 
the main cluster does not return to a state of thermal equi- 
librium at the end of the interaction since nothing prevents 
the central regions from cooling and the system moves to- 
ward extremely low values of luminosity and temperature. 
Fig. 1 191 illustrates the final structure of the remnants. In 
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Figure 20. X-ray surface brightness maps showing the late 
phases of the evolution of run l:6v3000. Two successive times 
are considered, when the centers of the projected total mass dis- 
tribution (white contours) are ~ 2500 kpc (top panel, Time = 2.4 
Gyr) and ~ 5500 kpc (bottom. Time = 4 Gyr) apart. Logarithmic 
color scaling is indicated by the key to the bottom of the figure, 
with violet corresponding to 0.9 keV and white to 86 keV. 



the case of the two highest velocities the final time is chosen 
in such a way that the distance between the centers of the 
two systems is about 1000 kpc larger than the sum of the 
virial radii at Time — 0. This occurs at Time = 2.2 and 4 
Gyr in the case of 1:6 and l:6v3000, respectively. The low 
velocity run l:6v2000 is analyzed at the time correspond- 
ing to the first apocenter when the core of the sub-cluster 
is close to the virial radius of the main cluster. The colli- 
sionless component of the main cluster (top left panel) in all 
cases is not substantially affected by the interaction while 
the sub-cluster dark matter halo appears to be strongly 
perturbed and retains the original spherical symmetry only 
within its scale radius. The sub-cluster central density profile 



(top right) is vertically shifted downward without a signifi- 
cant change of slope while for r > 0.2rvir the loss of material 
becomes more significant and in the case of low velocity en- 
counters {v = 2000 and 3000 km s~^) with mass ratio 1:6 the 
profile shows a large jump of more than one order of mag- 
nitude between O.Srvir and 0.5rvir- Beyond the scale radius 
the isodensity contours appear elongated and show a large 
plateau (see also Fig. [20]) associated to tidally stripped ma- 
terial. The evolution of gas density profiles is represented by 
the two panels on the bottom row of Fig. [19] In the 1:6 runs 
the interaction affects the central slope of the main cluster 
which becomes shallower while the 1:3 encounter produces 
the largest deviation from the initial values, with the final 
density profile shifted down by ^ 25%. As mentioned ear- 
lier, a not negligible part of the gas within the virial radius 
of the main cluster originally belonged to the sub-cluster 
and was subsequently stripped by ram-pressure during the 
central phases of the interaction. The density profiles of the 
stripped gas are drawn in blue for the different runs. In the 
case of l:6v2000 only the gas outside the virial radius of 
the sub-cluster is considered. Both l:6v3000 and l:6v2000 
are characterized by a large fraction of sub- cluster gas lost 
to the core of the main system, with a flat (l:6v2000) or 
even positive (l:6v3000) central slope, while the high veloc- 
ity run 1:6, despite the larger ram-pressure values, has less 
time to deposit gas in the central regions and shows a clear 
cut-off for r < O.lrvir- At larger radii, baryonic material is 
still accreting onto the remnant. Indeed, part of the main 
cluster gas has been pushed out by the bow shock and is 
now falling back into the cluster potential together with a 
fraction of the material lost by the sub-cluster. Curves in 
the bottom right panel of Fig. [19] represent the sub-cluster 
gas density profiles. The low velocity sub-clusters l:6v3000 
and l:6v2000 retain less than one tenth of their initial gas 
within O.lrvir while the density profiles in the external re- 
gions drop by a factor ^ 5. In general, the encounter flattens 
the gas density profile of the bullet core and this effect is ev- 
ident in the cooling simulation as well. The 1:6 bullet at 
Tiwie = 2.2 Gyr has lost almost all the gas within Vvir- Half 
of this baryonic material was stripped from its dark halo 
by ram-pressure during the central phases of the interaction 
and now accelerates towards its dark matter counterpart. It 
will be accreted by the sub-cluster halo at later times. 

Fig. [19] compares dark matter and gas density profiles 
at times when the main and sub-cluster are close to a state 
of virial equilibrium. Despite the changes observable in 
the gas density profiles and the tidal stripping affecting 
the bullet beyond its core radius, the central slope of the 
collisionless compone nt seems to be not stron gl y per turbed 
by the interaction (|Kazantzidis et al.1 l2004l . l2006h . The 
situation changes comparing the density profiles of the 
initial systems with those of the main and sub-cluster at 
the present time, when the bullet is located at almost 
one third of the host virial radius. As illustrated in the 
two upper panels of Fig. I21[ the density profile of the 
interacting dark matter halos increases in the inner regions 
and shows a decrement beyond OAvvir- This behavior is 
similar for different models and therefore independent of 
the orbital details, of the mass ratios and is seen for both, 
the main and the sub-cluster. The two-body relaxation time 
scale beyond the softening radius is longer than 10 "^ Gyr 
(JBovlan-Kolchin k MailiooJ : lArad k Johanssonll2005l ) . thus 
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implying that dark matter density profiles are not affected 
by numerical relaxation. As noticed in Section 3, fitting the 
projected masses around the centers of main and sub-cluster 
at the present time we would get more concentrated or more 
massive systems with respect to the real ones. The last 
two panels of Fig. [21] represent gas density profiles at the 
present time. The loss of gas in the central regions of the 
sub-cluster is expected as a result of ram-pressure stripping. 
But interestingly also the inner 200 kpc of the main cluster 
are completely devoid of gas originally belonging to the 
main cluster. Later on, after the displacement of the main 
cluster gaseous core, the baryonic material stripped from 
the bullet replenishes the central regions of the host system. 
As soon as the bullet moves toward the outskirts of the 
main system, the main cluster gas has time to collapse again 
into the center of the potential (bottom left panel of Fig. [19]). 



According to iHavashi fc White! (|2006l l the probability 
to find a cosmological configuration where the most massive 
sub-halo has a velocity larger than Vsub is well fitted by the 
function: 



logf{> Vsub) = -( 



Vsub /V200 \3.3 

1.55 ) ' 



(7) 



where ^7200 is the virial velocity of the main cluster. Adopt- 
ing for Vsub the present time velocity of the sub-cluster dark 
matter component (in Table [TJ we find a probability of 
5 X 10-^"^, 4.3 X 10-^ and 0.036 for runs 1:6, l:6v3000 and 
l:6v2000 respectively. 



5 CONCLUSIONS 

We used high resolution N-body/SPH simulations to per- 
form an extensive parameter study of the "bullet" cluster 
system 1E0657-56. Collisions of two NFW halos with hot 
isothermal gas components in hydrostatic equilibrium were 
studied, adopting initial relative velocities of 5000, 3000 and 
2000 km s~^, which is 3.8, 2.3 and 1.5 times the main clus- 
ter's virial velocity, respectively. We varied masses, mass 
ratios, impact parameter, baryonic fraction and concentra- 
tions. Most of the runs are adiabatic. Radiative cooling is 
included in two cases, adopting a standard cooling function 
for primordial gas. We analyzed the projected properties of 
the system at a time where the distance between the cen- 
ters of the mass distribution associated with the main and 
sub- cluster is comparable to the value provided by lensing 
observations. 

We have shown that: 

• Most of the main features in the observed X-ray maps 
are not well reproduced by encounters with zero impact pa- 
rameter. Indeed, depending on the relative velocity of the 
bullet, a perfectly head-on encounter either destroys the X- 
ray peak associated with the main cluster or does not pro- 
duce a significant displacement between gas and dark mat- 
ter in the main system. Moreover, in a pure head-on colli- 
sion, the gaseous core of the main cluster is displaced along 
the line which connects the centers of the two total mass 
distributions, contrary to observations. An impact param- 
eter corresponding to the core radius of the main cluster 
gas distribution provides enough ram-pressure to produce 
a displacement comparable to observations and introduces 



asymmetries in the main cluster emissivity map, similar to 
those detected in X-ray. 

• A low concentrated (c = 4) main cluster does not sur- 
vive the collision and its X-ray emissivity peak is destroyed. 

• Encounters with mass ratios as large as 1:3 do not 
match the observed X-ray morphology and the size of the 
projected temperature jump across the shock discontinu- 
ity. Introducing cooling in the simulations, the temperature 
peak is cooled to a value comparable with the observed one. 
On the other hand, lensing mass reconstructions seem to 
suggest even larger mass ratios. This remains an open ques- 
tion. Gravitational lensing analysis of the simulations pre- 
sented in this paper will be part of a follow-up paper (Mas- 
tropietro, Maccio' & Burkert in prep). 

• A pure cooling model neglecting energetic feedback still 
does not give a realistic description of the interacting sys- 
tem. What we are witnessing is overcooling in the central 
regions following a thermal instability in the early phases 
of the interaction. Nevertheless, cooling simulations provide 
a significant reduction of the temperature peak across the 
shock indicating that a more realistic treatment of the gas 
physics make models with higher mass ratios more realistic. 

• The choice of a different gas fraction does not affect 
significantly the results. 

• The morphology of X-ray maps is best simulated by adi- 
abatic 1:6 encounters. In particular the run l:6v3000, with 
initial velocity v = 5000 km s~^ (v ^ 2570 km s~^ in the 
center of mass system of reference) and present time dark 
matter velocity v ^ 3100 km s~^ (again in the center of mass 
rest-frame), reproduces most of the main X-ray features: the 
peculiar morphology of the X-ray emission associated with 
the main cluster, the relative surface brightness between the 
main and the sub-cluster, the shape of the bow shock and 
of the contact discontinuity. Decreasing the relative velocity 
(l:6v2000 has a present time dark matter velocity which is 
almost 300 km s~^ smaller than that used by Springel & 
Farrar 2007) the bullet becomes much less bright with re- 
spect to the center of the host system in disagreement with 
observations, although the displacement in the two clusters 
is better reproduced. The high velocity run 1:6 produces a 
contact discontinuity much broader than observed. 

• A significant fraction of the X-ray emission next to the 
center of the main cluster is associated with gas stripped 
from the external regions of the bullet. In the 1:6 runs the 
main cluster X-ray peak presents two distinct components: 
a compact strongly emitting region associated with the dis- 
placed core of the host system and a more diffuse component 
spatially coincident with gas stripped from the sub-cluster 
during the central phases of the interaction. The relative 
luminosity of this secondary component increases with de- 
creasing bullet velocities. 

• The observed difference in line of sight velocity is com- 
patible with an inclination of about 10 degrees of the plane 
of the encounter. A bottom-down inclination of our simula- 
tions keeps both the clusters close to the plane of the sky and 
does not produce essential changes in the the X-ray maps. 
But it affects the morphology of the main cluster peak and 
the luminosity of the secondary peak with respect to that 
associated with the main cluster core. 

• The projected temperature jump across the shock dis- 
continuity gives important indications about the nature of 
the encounter. Indeed, both, the height and the thickness 
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Figure 21. Dark matter (upper panels) and gas (bottom) density profiles of the main (left panels) and sub-cluster (right). Initial values 
and profiles at the present time are shown. Radius is scaled to the virial radius r^ir of the dark matter distribution. 



of the peak change, with the peak becoming broader for 
decreasing bullet velocities due to the lower pressure, exer- 
cised by the pre-shock gas after the central phases of the 
encounter. Assuming an adiabatic equation of state, both 
runs l:6v3000 and l:6v2000, match quite well the observed 
temperature jump, while for larger velocities (1:6) the peak 
becomes too narrow and pronounced compared to observa- 
tions. 

• Temperature maps reveal so me interesting feat u res. I n 
particular, as already found by ISpringel fc Farrarl (|2007l ), 
the bullet remains relatively cold despite of its dominant 
X-ray emission. For a short time, a high temperature re- 
gion appears next to the center of the main cluster mass 
distribution immediately after the central phases of the in- 
teraction. This feature - still visible at the present time - is 
not spatially coincident with the main cluster gaseous core. 
It is associated with hot gas stripped from the sub-cluster 
within the core of the host system. Its location partially cor- 



responds to the X-ray diffuse peak, while its eastern com- 
ponent co uld be related to the high temperature region ob- 
served by iMarkevitch et al.l (|2002l l southeast of the main- 
cluster peak. 

For a selected sub-sample of runs we performed a 
detailed three dimensional analysis following their past, 
present and future evolution for four Gyr. 

• Only the low velocity bullet l:6v2000 is actually bound 
to the main cluster at the end of the simulation while in 
l:6v3000 and 1:6 the sub-cluster has a velocity larger than 
the escape velocity from the host system. 

• The velocity of the gaseous component of the bullet 
starts diverging from that of its dark matter counterpart 
before the point of closest approach between the two clus- 
ters, whe n the gas is slowed down by ram-pressure. As pre- 
dicted bv lMilosavlievic et al.l (|2007f ), due to the drop in ram- 
pressure after the central phases of the interaction and grav- 
itational acceleration by its dark halo counterpart, at the 
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present time the gaseous bullet is not necessarily slower than 
the dark matter halo at the present time. We explored this 
question in greater details and found that the behavior of the 
sub- cluster gas after the point of closest approach strongly 
depends on the initial velocity of the bullet. For relatively 
low velocities (and ram-pressure values, run l:6v3000 and 
l:6v2000) the gaseous bullet is accelerated towards its dark 
matter counterpart as it leaves the core of the main cluster. 
At the present time the gaseous bullet moves as fast as its 
dark matter halo in l:6v3000 and even 500 km s~^ faster in 
l:6v2000. For larger encounter velocities (1:6) ram-pressure 
is more effective in slowing down the bullet even beyond the 
core radius of host cluster and the gaseous bullet is always 
slower than dark matter. 

• It ha s been already noticed by iMilosavlievic et al.l 
(12007*) and'Springel k FarraJ (|2007l ) that the sub-cluster ve- 
locities do not coincide with the shock velocity Vs as mea- 
sured by observers. In the case of the three runs (l:6v3000, 
l:6v2000 and l:6clfg) which best reproduce the projected 
temperature jump, the Mach number and shock velocity de- 
termined using the Rankine-Hugoniot conditions across the 
shock disco ntinuity show a good agreement with the values 
provided bv iMarkevitchI (|2006). The two dimensional shock 
velocity calculated by tracking the shock position as func- 
tion of time is consistent with these values. We find that, 
although after the point of closest approach Vs is always 
larger than the velocity of the sub-cluster dark matter halo, 
the difference depends on the model and becomes less signif- 
icant for higher velocities. In particular in the high velocity 
run 1:6 the present time shock velocity is only 6% higher 
than the dark matter bullet. 

• The collision produces a drastic increase in luminosity 
and temperature. The highest velocity impacts are associ- 
ated with the largest increase in temperature (and shocks) 
and the smallest peaks in luminosity due to the fact that in 
high velocity encounters a significant fraction of the baryons 
stripped from the bullet are deposit e d at l arge radii within 
the main cluster. iMarkevitch et al.l (|2002l l found that the 
bullet c luster lies on the Lx — T relation for local clusters 
(JMarkev itch 1998) at very high temperature. We followed 
the evolution of the main cluster in the Lx — T diagram 
during the different phases of the interaction and find that 
it moves roughly parallel to the Lx — T relation for nearby 
clusters. The maximum temperature is reached at the point 
of closest approach and the peak in luminosity immediately 
afterwards, with a small delay which produces a knot in the 
diagram. The branch of the curve associated with the post 
core-core interaction is still parallel to the Lx — T relation 
but shifted towards smaller luminosities with respect to the 
early increasing branch, due to the loss of gas beyond the 
virial radius of the main cluster. 

• After the encounter, as soon as the bullet is close or 
beyond the virial radius of its host system, the dark matter 
density profile of the main cluster does not deviate any- 
more significantly from the original one, while the gas pro- 
file becomes shallower in the central regions, with a large 
decrement in density observed only in run l:3clfg. The situ- 
ation changes drastically if we compare the density profiles 
of the initial systems with those at the present time, when 
the bullet is still located well within the virial radius of the 
main cluster. In this case the interacting systems are not 
in virial equilibrium and the dark matter densities of both. 



the main and sub-cluster, increase in the inner regions and 
show a decrement beyond 0.4 Vvir- Comparing these present 
time profiles with NFW halos would give wrong estimates of 
the halo parameters. At a time corresponding to the present 
configuration the center of the host system is completely de- 
void of main cluster gas as a result of the displacement of 
the main cluster gaseous core by the bullet. On the other 
hand, the baryons stripped from the bullet have replenished 
the central regions of the host. This is only a temporary 
situation however since a few Gyrs later the main cluster 
gas falls back and becomes again the dominant component 
in the central regions, although a significant fraction of gas 
stripped from the bullet is still present. The final sub-cluster 
dark matter density profiles seem significantly affected by 
the interaction beyond their scale radius, where the isoden- 
sity contours are elongated and show a plateau. 
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